Glassy behavior in a binary atomic mixture 
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We experimentally study one-dimensional, lattice-modulated Bose gases in the presence of an 
uncorrelated disorder potential formed by localized impurity atoms, and compare to the case of 
correlated quasi-disorder formed by an incommensurate lattice. While the effects of the two disorder 
realizations are comparable deeply in the strongly interacting regime, both showing signatures of 
Bose glass formation, we find a dramatic difference near the superfluid-to-insulator transition. In 
this transition region, we observe that random, uncorrelated disorder leads to a shift of the critical 
lattice depth for the breakdown of transport as opposed to the case of correlated quasi-disorder, 
where no such shift is seen. Our findings, which are consistent with recent predictions for interacting 
bosons in one dimension, illustrate the important role of correlations in disordered atomic systems. 

PACS numbers: 67.85.Hj ; 67.85. Fg ; 61.43.-j ; 03.75.Kk 



The presence of disorder is inherent to solid state sys- 
tems, and it has profound effects on transport in a variety 
of contexts, ranging from electron conductivity in metals 
to dirty superconductors [1]. Quantum gases in optical 
lattices [5] can, with a high degree of experimental con- 
trol, elucidate the role played by disorder in a number 
of physical phenomena. Recently, Anderson localization 
of matter-waves in disordered potentials has been ob- 
served for non-interacting gases [H |4j , and additionally, 
the reemergence of superfluidity due to repulsive inter- 
actions [5]. Discerning the (at times) competing roles 
of disorder and interactions is key to the understanding 
of Bose-glass behavior [5HS] in strongly interacting disor- 
dered systems. Ultracold atomic systems O [10] should 
provide a versatile testbed to aid in this endeavor [11] . 

Previous studies of ultracold atoms in disordered po- 
tential landscapes, generated by optical fields, have gen- 
erally suffered from strong correlations of the disorder 
that decay over length scales greater than either the 
healing length of the superfluid or the lattice spacing. 
This is true for both speckle potentials |31 [TU] that 
are diffraction-limited to structures on the order of the 
generating laser field's wavelength, and quasi-disordered 
bichromatic lattices [l] [S] |9] , which over large distances 
exhibit perfect correlations that may make them rather 
unsuitable for the realization of true disorder. To circum- 
vent these limitations, recently it has been proposed jl2|- 
I15j to use atomic impurities, which can be confined to 
regions much smaller than a lattice spacing, to act as 
point-like defects. 

The sudden quench of an atomic impurity field act- 
ing on mobile particles has been proposed [TSl [15] for 
the study of dynamical, out-of-equilibrium response to 
disorder. Alternatively, theoretical studies [TH [TB] have 
shown that even a slow "freeze-out" of the tunneling of 
one species from an initially homogenous mixture can 
lead to metastable "quantum emulsion" states, typified 
by local separation between frozen and mobile atoms 
(for repulsive interactions) and displaying properties sim- 



ilar to an equilibrium Bose glass. The study of quan- 
tum emulsions may help shed light on the coherence- 
loss mechanism in a number of experiments involving 
mass-imbalanced atomic mixtures, both for the boson- 
boson jTj and boson-fermion [151 112] cases. 

Here, we report on experimental studies of interact- 
ing one-dimensional (ID) Bose gases in the presence of 
disorder. We study the effects of uncorrelated disorder 
formed by atoms of an auxiliary spin state "frozen" to 
sites of an incommensurate lattice, and compare to the 
case of correlated quasi-disorder from an incommensurate 
bichromatic optical lattice. While both disorder types 
drive strongly interacting samples into an apparent Bose 
glass state, a large difference is seen for intermediate in- 
teractions, where we find that uncorrelated disorder has 
a dramatic effect in driving the system towards an in- 
sulating state. Our observation of enhanced localization 
for a more random disorder is consistent with recent the- 
oretical predictions for interacting bosons in ID [20] . 

The cartoon in Fig. [l] (a) qualitatively depicts our one- 
dimensional systems of lattice-trapped bosons with em- 
bedded impurities. After creating a homogeneous spin 
mixture, we slowly freeze the impurity atoms to a deep 
state-selective lattice of incommensurate spacing. For 
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FIG. 1. Disordered one-dimensional Bose gases, (a) To cre- 
ate a disordered impurity-field, half the atoms of a lattice- 
trapped Bose gas are converted to an auxiliary spin state 
(red) and localized to a state-selective incommensurate lat- 
tice (dashed green), (b) Alternatively, a weak incommen- 
surate lattice (green) is superimposed onto a lattice-trapped 
Bose gas. In both cases, the secondary potential (atomic or 
optical) causes site-dependent energy shifts Si and site-to-site 
energy differences A,. 
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FIG. 2. (a,b) Histograms of calculated A-distributions for the 
cases of the mixture (assuming either or 1 impurities per 
site) and an incommensurate lattice of depth s' = 1. Both 
are for a primary lattice depth s = 6 {U/Er = 0.4). (c,d) The 
autocorrelation function Xj = {AiAi+j)i/(AiAi)i of the A- 
distributions in (a,b), as a function of the site-to-site distance 
j. The averaging is over 1000 sites. 



comparison to a well-studied case of correlated disorder, 
we also study bosons in an incommensurate bichromatic 
lattice system [H IB H] , as depicted in Fig. [l] (b) . 

In both cases, the dynamics of the mobile atoms may 
approximately be described by the Bose-Hubbard Hamil- 



tonian (BHH) [71 [H] iJ = -tJ2^ia^, 
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T '^i ^i("-i~l)+Si '^jSji where t and U are the tunneling 
and interaction energies of the mobile atoms, and Oi, al, 
and fii — alcii are the annihilation, creation, and number 
operators for particles at lattice site i. The disorder- 
ing potentials will generally have two effects - to slightly 
modify the on-site wavefunctions of the atoms and to 
cause random site-dependent energy shifts e^. The first 
effect will lead to a spread of site-dependent values for 
both t and U (with a multi-band treatment necessary for 
strong perturbations). The second leads to random site- 
to-site energy differences A^ = Si — Ej+i, which define the 
resonance conditions for single-particle intersite tunnel- 
ing. In general, the many-body character will be defined 
by a competition between energy scales i, U, and A, with 
some dependence on the details of the A-distribution [7] . 
The two disorder potentials of Fig. [T] lead to quite dif- 
ferent A-distributions. In Fig. [2] we plot calculated his- 
tograms for (a) a fifty-percent impurity mixture and for 
(b) a weak incommensurate lattice (depth s' — 1, see 
experimental description below). While both distribu- 
tions are continuously filled and extend beyond ±A/f/, 
the details differ considerably. The impurity distribution 
is peaked about A/U = due to adjacent impurity-free 
sites, while the bichromatic lattice distribution is peaked 
at the outer-bounds of the distribution. In Fig. l2] (c,d) 
we plot the normalized autocorrelation function Xj — 
(AjA.i+j)i/(AiAj)j. The perfectly regular correlations 
in Fig. ^ (d) are a known attribute of quasi-disordered 
incommensurate lattices [Ml [55]. In contrast, for the 



atomic impurity field - combining the irregular spacing 
of bichromatic lattices and the irregular, probabilistic fill- 
ing of binary disorder - off-site correlations are strongly 
suppressed. This difference can have profound effects, 
as the localization properties of a system will generally 
depend both the strength and correlation length of the 
disorder potentials [53j. An extreme case can be found 
for non-interacting particles in ID, where random disor- 
der leads to Anderson localization for any finite disorder 
strength A 7^ 0, while incommensurate lattices induce 
localization only beyond a critical lattice depth ^l HI] ■ 

To briefly describe our experimental system, we be- 
gin as in [5S] with an optically trapped Bose-Einstein 
condensate of ^^Rb atoms. In 200 ms we load an ar- 
ray of isolated, one-dimensional tubes formed by the in- 
tersection of two optical lattices. These lattices are of 
period d = A/2 and depth 40 Er (with A = 1064 nm, 
Er — (/i/A)^/2m, Planck's constant h, m the atomic 
mass). The atoms are trapped along the tube axis z 
by a nearly harmonic potential of trapping frequency 
W2/27r = 80 Hz. A lattice along z, also of period d and 
with variable depth s (in units of the recoil energy Er)^ is 
smoothly ramped up within 100 ms. This primary lattice 
serves to define the sites (index i) and parameter values 
(t, U) of the BHH. 

Initially, the tubes contain only atoms in the |-F, mp) = 
|2, — 2) hyperfine ground state. To create atomic impuri- 
ties as in Fig. IT] (a), a fraction (/imp) of the total popula- 
tion of 8 X 10'* atoms is transferred to the 1 1, — 1) state via 
a microwave Landau-Zener sweep. The impurity atoms 
are loaded into a completely state-selective lattice [25] 
along z in 20 ms. This lattice has spacing d' — A'/2, with 
A' = 785 nm, and the impurity atoms are deeply local- 
ized at 20 E'j^ (recoil energy E'j^ = {h/\'f/2m - 1.8£;j^.). 
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FIG. 3. Disappearance of excitation gap due to disorder, (a) 
Visibility as a function of amplitude modulation frequency 
(normalized to U/Er — 0.53 for s = 14), in the absence of 
disorder for s = 9 and s = 14 (open and filled black circles). 
Fit lines are two Gaussians on a linear slope, (b) For s = 
14, with atomic impurities, open purple squares and filled 
red squares represent /imp — 0.1 and 0.5, respectively, (c) 
For s — 14, with an incommensurate lattice of depth s' — 1 
(orange triangles) and no impurities. Solid error bars are 
statistical over several runs, while dashed are estimated errors 
for individual runs (120% of maximum statistical error). 



We typically study an equal mixture (/imp = 0.5) of 
frozen and mobile atoms, which we shall refer to as "the 
mixture". For the alternate disorder implementation of 
Fig. [l] (b), we begin with a sample of 4 x 10'* atoms in 
the state |2, —2) (i.e. same number of mobile atoms). We 
then ramp up a secondary lattice in 20 ms, of spacing d! 
and variable depth s' x E'^^, onto the |2, —2) atoms. Due 
to the external trapping potential, our system is not ho- 
mogeneous, but can be characterized by a typical central 
filling factor of n ~ 3 (total) atoms per site. 

We begin our investigation into the effects of disor- 
der by measuring excitation spectra, which relate most 
directly to the distribution of site-to-site energy shifts. 
Along with a finite compressibility [35], a gapless exci- 
tation spectrum is a characteristic feature distinguishing 
a disordered Bose glass state from a (homogenous) Mott 
insulator. We measure the excitation spectra by per- 
forming amplitude- modulation spectroscopy [9l |27l |29] of 
the primary z-lattice at driving frequencies w„iod/27r [50] . 
For the disorder- free case (A^ « 0) in Fig. [s] (a), with 
the sample chosen to be deep into the ID Mott regime 
(s = 14 , U /t w 66), the excitation spectrum exhibits res- 
onant structure. The resonance positions are consistent 
with the excitation of particle-hole pairs at U /h [571 [SI] 
(and 2U /h due to either higher-order processes or excita- 
tion at the edge of Mott domains [171I3I]). In contrast, 
for both the atomic impurity mixture (Fig. [31(b)) and for 
an incommensurate lattice of depth s' = 1 (Fig. [3] (c)) 
having comparable A-distribution bounds, we observe 
flat excitation spectra (cf. [3]). These observations are ex- 
pected [7] for broadly-filled A-distributions with bounds 
Aniax > U , and are consistent with the system being 
in a Bose glass state (future compressibility measure- 
ments [26j should allow for the disambiguation between 
a true Bose glass [311] and a disordered Mott state ^S]). 

While the observed spectral properties are consistent 
with Bose glass formation, transport measurements are 
necessary to confirm insulating behavior. Here, we study 
the effects of disorder in the transition region between 
superfluid and insulator, determining the critical lattice 
depth at which the systems become insulating through 
the study of localization and transport. In regard to the 
former, the momentum-peak width of a released sam- 
ple (related to the inverse correlation length £^~^ of the 
sample in-situ [35]), exhibits a sudden increase accom- 
panying an abating superfluid fraction and loss of off- 
diagonal long-range order |33j. As for the latter, it has 
been shown [27l [34] that the response to an applied im- 
pulse dies away upon entering the strongly correlated 
regime, and can serve as a signature of insulating be- 
havior [Tni[2ni- 

To study transport, we look at how the system evolves 
after an applied impulse. A magnetic-field gradient along 
z is pulsed on for duration of T = 1.2 ms and applies a 
variable force F ranging from to Fmax/m = 1.2 m/s^, 
resulting in an impulse I = F x T. As illustrated in 
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FIG. 4. (a) The response to impulse is determined by a 
straight-line fit (with slope a, normalized to the case of free 
atoms) to the dependence of velocity on applied impulse / 
(profiles shown for disorder-free and impurity mixture cases; 
recoil velocity vr = h/2md). (b) a versus lattice depth s 
for bosons without disorder (black circles) and with atomic 
impurities (/imp = 0.5, red squares). Lines and surrounding 
shaded regions are fits to the data with confidence regions 
(1 s.d.) of a linearly decaying response, with no response 
(a = 0) beyond a depth Sc- (c) Similarly, but for an incom- 
mensurate lattice (s' = 3, green diamonds), with disorder-free 
data reproduced for comparison. The inset shows values of Sc 
versus s' as determined by fit intercepts, (d) The momentum- 
peak width a {l/\/e half-width) is determined by a symmetric 
multi-Gaussian fit to time-of-flight interference patterns. The 
shaded regions about v/vji = 0(±1) are used to count A''o(±i) 
for visibility measurements of Fig. [s] (e) a versus s for ID 
bosons without disorder and for the mixture (colors/symbols 
as in (b)). The straight line for the disorder-free data is a 
fit of a linear increase beyond a depth Sc- The first few data 
points for the case of atomic impurities are connected as a 
guide to the eye. (f) Similarly, but for an incommensurate 
lattice with s' = 2 (blue triangles) and with disorder-free 
data reproduced. 



Fig. H (a) , we characterize the response as a function of 
/ (with slope a) by monitoring the center-of-mass ve- 
locity along z in time-of-flight absorption images follow- 
ing a brief (~1 ms) ramp-off of the z-lattice. We access 
the momentum-peak width as in [SJ [57] by releasing the 
atoms in time-of-flight following a 50 /iS lattice ramp up 
to s = 20 and a gravitational phase-shift along z (with- 
out impulse or lattice ramp-off). We then determine the 
peak width cr by a fit to the profile of symmetric diffrac- 
tion peaks on top of an incoherent background, as shown 
inFig.|4](d). 

We find that for the impurity mixture, the mobile 



atoms are more easily driven towards insulating behavior. 
Similar to [29], we determine the critical point at which 
the atoms become unresponsive to impulse by fitting a 
linear decay to the response a as a function of lattice 
depth. As shown in Fig. l4] (b) , this fit to the transport 
measurements yields a critical depth Sc = 3.0±0.5. While 
we do not fit the peak width data, a kink near Sc ~ 2 can 
be observed in Fig. H (e) . These values are roughly half 
of those measured for the disorder-free case, with values 
of Sc = 6.2 ± 0.5 and 5.7 ± 0.5 based on impulse response 
and peak width, respectively. This value of Sc « 6 is 
close to our expectations based on the mean Lieb-Liniger 
parameter [35] (7 = 0.6) of our ID gases prior to lattice- 
loading, with an estimate based on the Bose-Hubbard 
(Sine-Gordon) model predicting s^ = 6.2 (4.7) [^ . 

For the incommensurate lattice, transport data for s' = 
3 is shown in Fig. l4] (c) . In this case no shift of the critical 
depth is seen, and as shown in the inset, this is the case 
for all incommensurate lattice depths considered (s' < 3, 
where we restrict to s' /s < 1 in the transition region to 
maintain the perturbative nature of the incommensurate 
lattice). The momentum-peak width data mirrors this 
lack of a shift of the transition point (inset: for all depths 
considered s' < 2). 

While both incarnations of disorder resulted in an ap- 
parent Bose glass state for very deep lattices, a clear dif- 
ference was seen in their effect on more weakly interacting 
samples. In attempting to account for the observed dif- 
ference, a natural consideration is the disparity in their 
correlations [xj, cf. Fig.[2](c,d)]. In general, one expects 
that the less correlated the disorder, the more enhanced is 
the localization [53]. For interacting bosons in ID, it has 
been shown theoretically [201 that the localization tran- 
sition occurs for a more weakly interacting gas (larger 
values of the Luttinger exponent K or lower values of 
the Lieb-Liniger parameter 7 [24] ) in uncorrelated disor- 
der than for correlated disorder. Our observations of a 
sizeable shift of the transition point for impurities and a 
negligible shift for an incommensurate lattice are thus in 
qualitative agreement with expectations based on their 
dissimilar correlation properties. 

Also relevant to our observations is the reduction of 
phase-space density in the presence of localized impuri- 
ties, as well as of the atomic density due to the dynamical 
formation of impurities from the mobile species. The first 
effect has been shown [36] to be responsible for adiabatic 
heating and loss of coherence in recent Bose-Fermi mix- 
ture experiments [TO], due to reduced entropy following 
a reduction in effectively occupiable sites, both for at- 
tractive and repulsive interactions. The second effect is 
more particular to "quantum emulsion" [TTj [HI US] ex- 
periments. Here, reduced density leads to more strongly 
correlated many-body states in ID and thus favors in- 
creased localization and insulating behavior. 

In conclusion, we have observed signatures of Bose 
glass formation in ID Bose gases with superimposed 



disorder, both for atomic impurities and for quasi- 
disordered bichromatic lattices. The two disorder types 
have dramatically different effects in the transition region 
between superfluid and insulator, with atomic impurity 
disorder inducing localization in much more weakly in- 
teracting gases. Our observation that a more weakly 
correlated disorder leads to enhanced localization is in 
qualitative agreement with recent theoretical predictions 
for interacting ID boson systems. 
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